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A new non-perturbative approach to quantum field theory — D-theory — is proposed, in which continuous 
classical fields are replaced by discrete quantized variables which undergo dimensional reduction. The 2-d clas- 
sical 0(3) model emerges from the (2 + l)-d quantum Heisenberg model formulated in terms of quantum spins. 
Dimensional reduction is demonstrated explicitly by simulating correlation lengths up to 350,000 lattice spacings 
using a loop cluster algorithm. In the framework of D-theory, gauge theories are formulated in terms of quantum 
links — the gauge analogs of quantum spins. Quantum links are parallel transporter matrices whose elements 
are non-commuting operators. They can be expressed as bilinears of anticommuting fermion constituents. In 
quantum link models dimensional reduction to four dimensions occurs, due to the presence of a 5-d Coulomb 
phase, whose existence is confirmed by detailed simulations using standard lattice gauge theory. Using Shamir's 
variant of Kaplan's fermion proposal, in quantum link QCD quarks appear as edge states of a 5-d slab. This 
naturally protects their chiral symmetries without fine-tuning. The first efficient cluster algorithm for a gauge 
theory with a continuous gauge group is formulated for the U(l) quantum link model. Improved estimators for 
Wilson loops are constructed, and dimensional reduction to ordinary lattice QED is verified numerically. 



1. Introduction 

The conventional approach to field quantiza- 
tion is to perform a path integral over classical 
fields. Wilson's formulation of lattice field theory 
maintains close contact with perturbation theory 
by following this approach. This has obvious ad- 
vantages, for example, in QCD when one wants 
to make contact with continuum perturbative re- 
sults. On the other hand, the close relation to 
perturbative methods may imply that it is un- 
necessarily difficult to probe the non-perturbative 
regime. Here we propose an alternative method 
to quantize field theories. Instead of using contin- 
uous classical fields, we work with discrete quan- 
tized variables. In scalar field theories these vari- 
ables are generalized quantum spins, and in gauge 
theories they are quantum links. It may seem 
paradoxical that models with continuous symme- 
tries can be expressed using discrete variables. 



Note, however, that — although a spin 1/2 has 
only two discrete states — the continuous 50(3) 
symmetry is represented exactly. D-theory gen- 
eralizes this structure to other symmetries. 

Why should this quantization procedure be 
equivalent to performing a path integral over clas- 
sical fields? A quantum spin operates in a finite 
Hilbert space. How can the full Hilbert space 
of a quantum field theory be recovered when the 
classical fields are replaced by analogs of quan- 
tum spins? As we will see, this requires specific 
dynamics: the dimensional reduction of discrete 
variables. This dynamics arises in quantum anti- 
ferromagnets as well as in quantum link models 
||, which are the gauge analogs of quantum spin 
systems [^],[| . The dimensional reduction of dis- 
crete variables leads to a new non-perturbative 
D-theory formulation of quantum field theory [Q. 
Quantized variables can be used e.g. to formulate 
0(N), SU(N)®SU(N) and CP(N) spin models, 
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as well as Abelian U(l) and non-Abelian SO(N), 
SU(NY and U(N) gauge theories, in particular, 
QCD |. 

D-theory is closely related to Wilson's lattice 
field theory. The classical fields of Wilson's for- 
mulation arise as collective excitations of a large 
number of discrete D-theory variables. D-theory 
is defined on a much finer lattice than the result- 
ing effective theory, and hence provides a micro- 
scopic structure underlying Wilson's theory. Our 
hope is that the microscopic structure simplifies 
non-perturbative calculations, and hence may be 
helpful to understand the dynamics. For exam- 
ple, due to the discrete nature of its degrees of 
freedom, D-theory can be simulated with very ef- 
ficient cluster algorithms. Also, the theory can 
be completely fermionized, which leads to an in- 
teresting approach to the large N limit. D-theory 
shares several interesting properties with the non- 
perturbative formulation of string theory, and can 
hence also be viewed as a daughter of M-theory. 

This paper is based on five talks. In section 
2, dimensional reduction of discrete variables is 
illustrated in the context of quantum antiferro- 
magnets. Section 3 describes a finite-size scaling 
analysis that confirms dimensional reduction of 
the (2 + l)-d quantum Heisenberg model to the 
2-d classical 0(3) model using a loop cluster algo- 
rithm. In section 4, (4+l)-d quantum link models 
are constructed and their dimensional reduction 
to ordinary 4-d Wilsonian lattice gauge theories is 
discussed. Section 5 describes a numerical study 
of 5-d non-Abelian gauge theories, confirming the 
presence of a Coulomb phase, which is essential 
for dimensional reduction. In section 6, the D- 
theory formulation of full QCD is obtained by in- 
cluding quarks as edge states on a 5-d slab. The 
first efficient cluster algorithm for a gauge theory 
with a continuous gauge group is constructed in 
section 7 for the U(l) quantum link model. Fi- 
nally, section 8 contains our conclusions. 

2. The 2-d 0(3) Model from Dimensional 
Reduction of Quantum Spins 1 

Quantum antiferromagnets like La2Cu04 and 
S^CuC^Ck) are the precursor insulators of high- 

1 Based on a talk given by U.-J. Wiese 



temperature superconductors. These materi- 
als are well described by the 2-d antiferromag- 
netic square-lattice spin 1/2 quantum Heisenberg 
model, which is formulated in terms of quantum 
spin operators S x (Pauli matrices) with the usual 
commutation relations 

[Si,Si\=iS xy e ijk Sl (1) 

The corresponding Hamilton operator 

H = J ^ S x ■ S x+ p,, (2) 

is defined on a 2-d square-lattice. Here we 
consider antiferromagnets, i.e. J > 0. The 
problem has a global SO(3) symmetry, because 
[H, J2 X = 0> i- c -i the Hamilton operator com- 
mutes with the total spin. The quantum statisti- 
cal partition function takes the form 

Z = Trexp(-/3H). (3) 

When Z is expressed as a path integral, the in- 
verse temperature (3 occurs as the extent of the 
Euclidean time direction. In the zero tempera- 
ture limit the quantum system resembles a 3-d 
classical model at infinite volume. The 2-d quan- 
tum Heisenberg model (^-|| as well as the exper- 
imental systems M exhibit long-range antiferro- 
magnetic order. Consequently, the 50(3) sym- 
metry of the corresponding 3-d classical model is 
spontaneously broken to SO(2). Then, due to 
Goldstone's theorem, two massless bosons — in 
this case two antiferromagnetic magnons (or spin- 
waves) — arise. Using chiral perturbation the- 
ory, the low-energy dynamics can be described 
by effective fields in the coset SO(3)/SO(2) = S 2 
|l(J. Hence, the magnon field is a unit vector s — 
the same field that appears in the classical 0(3) 
model. This is how the usual path integral vari- 
ables arise from D-theory. Due to spontaneous 
symmetry breaking, the collective excitations of 
many discrete quantum spin variables form an ef- 
fective continuous classical field. 

To lowest order in chiral perturbation the- 
ory, the effective action describing the low-energy 
magnon dynamics is given by 
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Here c and p s are the spin-wave velocity and the 
spin stiffness. The index /z extends over the spa- 
tial directions 1 and 2 only. At finite tempera- 
ture T, the 2-d quantum system resembles a 3- 
d classical model with finite extent (3 = l/T in 
the Euclidean time direction. At zero tempera- 
ture spontaneous symmetry breaking implies the 
existence of massless Goldstone bosons with an 
infinite correlation length £. At non-zero tem- 
perature the extent (3 of the extra dimension 
becomes finite. Then the system appears di- 
mensionally reduced to two dimensions, because 
(3 <C £. However, in two dimensions the Mermin- 
Wagner-Coleman theorem prevents the existence 
of interacting massless Goldstone bosons Jll| . As 
a consequence, in the 2-d 0(3) model a mass gap 
is generated non-perturbatively. 

Is the now finite correlation length large enough 
for dimensional reduction still to occur? This 
question was first answered by Chakravarty, 
Halperin and Nelson |{l2| , and then studied in 
more detail by Hasenfratz and Niedermayer p3[ . 
They used a block spin renormalization group 
transformation that maps the 3-d chiral pertur- 
bation theory model with finite extent f3 to a 2-d 
lattice 0(3) model. The 3-d magnon field is aver- 
aged over volumes of size (3 in the Euclidean time 
direction and (3c in the two spatial directions. 
Due to the large correlation length, the field is 
essentially constant over these blocks. The aver- 
aged field is defined at the block centers, which 
form a 2-d lattice, whose spacing (3c is much 
coarser than the microscopic lattice spacing a of 
the original quantum antiferromagnet. The effec- 
tive action of the averaged field resembles a 2-d 
lattice 0(3) model formulated in Wilson's frame- 
work. However, since the effective coarse lattice 
action is obtained by exact blocking of the origi- 
nal fine lattice degrees of freedom, its lattice ar- 
tifacts are entirely due to the microscopic lattice. 
This is important for the practical applicability of 
D-thcory. Using the 3-loop /^-function of the 2-d 
0(3) model together with its exact mass-gap [jl4| , 
Hasenfratz and Niedermayer JT3J generalized an 
earlier result of Chakravarty, Halperin and Nel- 



son |12(] to what we call the CH2N2-formula 



167T/9; 



exp(27r/3p s )[l 



Att(3 Ps 



Here e is the base of the natural logarithm. Since 
in the zero-temperature limit £ is infinitely larger 
than (3, the system undergoes dimensional reduc- 
tion to two dimensions when the extent of the 
extra dimension becomes large. The exponen- 
tial increase of £ in the zero-temperature limit is 
a consequence of asymptotic freedom of the 2-d 
classical 0(3) model. The factor 2ir in the expo- 
nent is the 1-loop coefficient of the corresponding 
/3-function, and 
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is the inverse coupling constant of the dimen- 
sionally reduced lattice model. Hence, in the 
low-temperature limit, the 2-d antiferromagnetic 
quantum Heisenberg model provides a new non- 
perturbative regularization of the 2-d 0(3) model. 
This D-theory formulation is entirely discrete, 
and still has the continuous SO(3) symmetry. 

The dimensional reduction of a (d + 1)- 
dimensional D-theory (in this case the (2 + 1)- 
d quantum Heisenberg model) to an effective d- 
dimensional Wilsonian lattice theory with lattice 
spacing (3c is illustrated in fig.l. The continuum 
limit of the effective lattice theory is reached as 
g = l/(3p s — > 0, and hence as the extent (3 of the 
extra dimension becomes large. Still, in physical 
units of £, the extent (3 <C £ becomes negligible 
in this limit. In the continuum limit, the lattice 
spacing (3c of the effective 2-d Wilsonian lattice 
0(3) model becomes large in units of the micro- 
scopic lattice spacing of the quantum spin sys- 
tem. Hence, D-theory introduces a discrete sub- 
structure underlying Wilson's lattice theory. Due 
to exact blocking the lattice artifacts are entirely 
due to the microscopic D-theory lattice. 

3. Continuous-Time Loop Cluster Simula- 
tion of the Heisenberg Model at Very 
Large Correlation Lengths 2 



2 Based on a talk given by B. B. Beard 
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d-dimensional ordinary lattice field theory 




Figure 1. Dimensional reduction of a D-theory: 
Averaging the (d + l)-dimensional effective field 
of the D-theory over blocks of size (3 in the extra 
dimension and (3c in the physical directions re- 
sults in an effective d-dimensional Wilsonian lat- 
tice field theory with lattice spacing (3c. 

Soon after the discovery of high-temperature 
superconductivity in doped lamellar copper ox- 
ides it was found that the undoped com- 
pounds are 2-d spin- 1/2 quantum antiferromag- 
nets. Through experimental, theoretical, and nu- 
merical efforts much progress has been made in 
the understanding of these systems. In partic- 
ular, detailed neutron scattering measurements 
H of the correlation length £ in the magnet 
Sr 2 Cu0 2 Cl 2 were found to be in apparent good 
agreement with the CH 2 N 2 -formula eq.@ — i.e. 
with 3-loop asymptotic scaling of the effective 2- 
d lattice 0(3) model resulting from dimensional 
reduction. However, neutron scattering mea- 
surements on higher-spin systems |g|Jl5| reveal 
a striking discrepancy with the CH 2 N 2 -formula. 
The puzzle concerning the applicability of 3-loop 
asymptotic scaling has recently been resolved by 
simulating the Heisenberg model at very large 
correlation lengths £ w 350,000a |lq| . 

In the context of relativistic quantum field the- 
ory, where the lattice spacing serves as an ul- 
traviolet cut-off that is ultimately removed, the 
question of asymptotic scaling is unphysical, be- 



cause it involves the bare coupling constant. At 
low temperatures a 2-d quantum antiferromagnet 
induces a lattice action for an effective 2-d 0(3) 
model with lattice spacing a' = (3c, which is much 
larger than the microscopic lattice spacing a of 
the quantum antiferromagnet, and which is de- 
termined by the physical temperature T = 1/(3. 
Hence, the question of asymptotic scaling be- 
comes a physical issue. A priori, it is unclear for 
what values of £ one should expect asymptotic 
scaling for the effective 2-d lattice action. How- 
ever, it is known that asymptotic scaling often 
sets in only at very large correlation lengths. For 
example, in the 2-d classical 0(3) model with the 
standard lattice action, 3-loop asymptotic scal- 
ing sets in at about £ « 10 5 lattice spacings Q . 
Hence, one might expect that the CH 2 N 2 -formula 
works only above £ w 10 5 a' = 10 5 (3c. How- 
ever, the induced effective lattice action is not 
the standard action. Indeed, the effective action 
results from exact blocking of the chiral pertur- 
bation theory action of the magnons. This action 
receives cut-off effects only from the microscopic 
lattice with spacing a. In the limit a — > the 
induced coarse lattice action would be a perfect 
action. This means that in D-theory the cut-off 
effects are due to the fine lattice spacing a, not 
due to the lattice spacing a' = (3c of the induced 
d-dimensional effective Wilsonian theory. Based 
on this argument, one expects that in the quan- 
tum Heisenberg model 3-loop asymptotic scaling 
sets in at £ « 10 5 a. Indeed, this is observed in 
the numerical simulations. 

How can one investigate correlation lengths in 
the range £ w 10 5 a? In the classical 2-d 0(3) 
model this was possible using finite-size scaling 
methods [[j^jjlT]]. To investigate the correlation 
length in the quantum Heisenberg model we use 
the same technique. The key observation of finite- 
size scaling is that the finite-volume correlation 
length £(L) of a periodic system of spatial size 
L x L is a universal function of £/L, where £ is 
the correlation length of the infinite system. In- 
verting this relation one writes 

£ = L/(l/z) = L/(£(L)/L). (7) 

Here z = L/£(L) is Fisher's finite-size scaling 
variable, which is a renormalization group invari- 
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ant measure of the physical volume. For the clas- 
sical 2-d 0(3) model the universal function / has 
already been determined very precisely ||il||l7f . 
Since at low temperatures the quantum Hcisen- 
berg model reduces to a classical 2-d lattice 0(3) 
model, we can use the same universal function to 
deduce infinite volume results from finite-volume 
correlation length data. It is important that 
eq.(^) assumes universal behavior, i.e. scaling, 
but not asymptotic scaling. Hence, we can test 
the CH 2 N 2 -formula without bias. 

The finite-size scaling function / is very sensi- 
tive to small changes in the finite- volume correla- 
tion length £(L). A small error in generates 
large uncertainties in the infinite-volume corre- 
lation length £. Hence, one needs a very accu- 
rate numerical method to determine £{L). For- 
tunately, for the quantum Heisenberg model a 
very efficient loop cluster algorithm exists p9|,[7[, 
which practically eliminates auto-correlations in 
the Monte Carlo data. Also, it is possible to con- 
struct improved estimators which drastically re- 
duce statistical errors. The only remaining sys- 
tematic error is due to the Suzuki- Trotter dis- 
cretization of Euclidean time. Recently, it has 
been realized that the discretization of Euclidean 
time is not necessary for path integrals of dis- 
crete quantum systems |^] . This completely elim- 
inates the systematic discretization error of previ- 
ous methods. Also, it greatly reduces storage and 
computer time requirements, and thus allows us 
to work at very low temperatures. 

We have performed simulations at inverse tem- 
peratures J/T — 0.5, 1.0, 3.5 measuring the 
correlation length directly in a large volume L « 
6£. For J/T = 4.0, 4.5, 12 we have used the 
finite-size scaling technique. The continuous time 
loop cluster algorithm was used in a single cluster 
version, and an improved estimator has been im- 
plemented for the staggered correlation function. 
Correlation lengths are extracted using the sec- 
ond moment method |lqjl7| . In all cases at least 
10 5 measurements have been performed. The nu- 
merical data for the infinite-volume correlation 
length £ are compared with experimental data 
and the CH2N2-result in fig. 2. The numerically 
accessible correlation lengths are more than three 
orders of magnitude larger than the experimental 
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Figure 2. The correlation length £/a: Experi- 
mental data (crosses) and numerical data (circles) 
are in apparent good agreement with the CH2N2- 
formula eq.(|5|) (dashed line). 

ones jj). 

Focusing on small effects invisible in fig. 2, fig. 3 
shows the deviation from 2-loop asymptotic scal- 
ing as a function of temperature. The correla- 
tion length data have been fitted simultaneously 
with previously obtained data for the staggered 
and uniform susceptibilities in cubic and cylin- 
drical HI space-time geometries. Indeed, we find 
that asymptotic scaling at the 3-loop level of the 
CH2N2-formula sets in at correlation lengths of 
£ » 10 5 a, while a 4-loop fit works well already at 
£ w 10 2 a. When the massgap is included as a fit 
parameter, its fitted value agrees with the exact 
massgap at the 2 percent level. In the light of our 
numerical data, the apparent good agreement of 
the CH2N2-formula with the less accurate spin 
1/2 experimental data at higher temperatures 
must be considered accidental. The large discrep- 
ancy between the CH2N2-formula and experimen- 
tal data for higher-spin systems ||[l5| probably 
arises, because again 3-loop asymptotic scaling 
sets in only at very small temperatures, which 
are inaccessible to experiments. An investigation 
of the spin 1 case using the continuous-time loop 
cluster algorithm is presently in progress. 

The computational effort for simulating the 
quantum Heisenberg model is compatible to us- 
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Figure 3. Analysis of the correlation length 
£ as the asymptotic scaling regime is ap- 
proached: £ is divided by the 2-loop result £2 = 
(ec/167r/5 s ) exp(27r/9 s /T), using the best-fit value 
p s = 0.1800(5). The 3-loop result (dashed line) 
is accurate only for £ « 10 5 a, while the 4-loop 
regime (solid line) begins at £ s=s 10 2 a. 

ing the Wolff cluster algorithm p3] directly in 
the 2-d 0(3) model. For gauge theories, on the 
other hand, no efficient cluster algorithm has been 
found in Wilson's formulation of the problem. 
Recently, the first efficient cluster algorithm for 
a gauge theory with a continuous gauge group 
has been constructed for the U(l) quantum link 
model 0|. If efficient cluster algorithms can be 
constructed also for non-Abelian quantum link 
models, numerical simulations of QCD will be- 
come more accurate than the ones using Wilson's 
method. 

4. QCD as a Quantum Link Model 3 

In Wilson's formulation of lattice gauge theory 
the parallel transporters u x ^ are classical SU(N) 
matrices defined on the links (x, //) of a 4-d hyper- 
cubic lattice. The standard Wilson action takes 
the form 

%] = - TT \- u x^ u x+^ u i+6^ ut xJ- ( 8 ) 

3 Based on the talks given by R. Brower and U.-J. Wiese 



By construction, this action is invariant under 
gauge transformations 

u xfl = exp(ia x • \)u x ^ exp(—ia x +p, ■ A). (9) 

In D-theory the Wilson action is replaced by a 
quantum link model Hamilton operator 

H = J ^P x ,„U x+ ^Ut + ^Ul v \. (10) 

Wilson's classical link matrices are turned into 
quantum link operators U x ^. They still are 
N x N matrices, but now their elements are non- 
commuting operators acting in a finite Hilbert 
space. Of course, we want to maintain the gauge 
symmetry of the problem. In quantum link mod- 
els, gauge invariance means that H commutes 
with local generators G x of gauge transformations 
at the site x, which obey the usual algebra 

[G a x ,G b y ] = 2i5 xy f abc G x . (II) 

Under a unitary transformation JT exp(ia x - G x ), 
representing a gauge transformation in Hilbert 
space, a quantum link variable transforms as 

U x , IJ , = Y]. cx P(~ i &y ' Gy)U x ^^eyij>{ia z ■ G z ) 

V z 

= exp(ia x • X)U Xtlx ex.p(-id x+(i ■ A). (12) 
This relation implies 

[G X , Uy^] — 5 X ,y+p,Uy tlJ/ )l — S x .y\Uy,^, (13) 

which can be satisfied when we introduce 

G x = ^ ^ (Rx-p.,ii + L x ^). (14) 

Here L x ^ and R x ,^ are generators of left and right 
gauge transformations of the link variable U Xtfl . 
They are 2(N 2 - 1) generators of an SU(N) L ® 
SU(N)r algebra on each link. The commutation 
relations of eq.([l^) imply 

[-Lx : fj,i Uy^jj\ — 8 X ^y8fiiyXU x ^^^ 

[Rx,lll Uy >V \ = 5 X .y5^ V U X ,^\. (15) 

In D-theory, the real and imaginary parts of the 
elements of the link matrices are represented by 



2N 2 Hermitean operators. Together with the op- 
erators L x ,ii and Rx,^, these are 4N 2 — 2 genera- 
tors. The above commutation relations are those 
of an SU (2N) algebra, which contains SU (N) l ® 
SU(N) R as a sub-algebra. Of course, SU(2N) 
has AN 2 — 1 generators, and indeed there is one 
more generator T x ^, which obeys 



Consequently, 



(16) 



(17) 



generates an additional £7(1) gauge transforma- 
tion 

u 'x,n = H exp{-iayG y )U XjfJ . e-np(ia z G z ) 
y * 
= cx_p(ia x )U x ^exp(-ia x+t _ t ). (18) 

Indeed, the Hamilton operator of eq.(|lO|) is also 
invariant under the extra U(l) gauge transforma- 
tions and thus describes a U(N) lattice gauge the- 
ory. The symmetry can be reduced to SU(N) by 
adding the real part of the determinant of each 
link matrix to the Hamilton operator 

H = J ]T Tr[U x ,,U x+ ^Ul + ^UU 
+ J'J2 [detU^ + detJZy. 

X,fJ, 



(19) 



Due to the discrete nature of the quantum link 
variables, the above algebraic structure can be 
expressed entirely in terms of anticommuting op- 
erators 

Tjij _Jt J 

f — X 



Rx,fj, 



T x ,fj, — /~^,( C x,+u C x,+U C X+fl. 



- m 4+a,-m)- ( 2 °) 



Here c l J ±„ and cj. ±„ are creation and annihila- 
tion operators of colored fermions living on a link 
that emanates from the site x in the indirection. 
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x,H 



Figure 4. QCD dynamics as a rishon abacus: The 
trace part of the Hamiltonian induces hopping of 
rishons of various colors around a plaquette. The 
determinant part shifts a color-neutral combina- 
tion of N rishons from one end of a link to the 
other. 

In quantum link models these fermions can be 
viewed as constituents of the gluons, which we 
call rishons. The rishon operators obey canonical 
anticommutation relations 

{ c x,±fii Cy,±t/} = &xyfi±ix,±v&iji 

{ c x,±tii ^y,±v) — { c x,±/ii ^y,±v) = 0- (21) 

It is easy to show that the rishon number 



£^\ u x+fi,-iJ, c x+p,-iJ + c x' +^ c k,+p)j (22) 



is conserved separately for each link. In rishon 
representation the determinant of a quantum link 
operator takes the form 

det U x ,„ = An4 )+ ^4j_ M ..< +M c£l M . (23) 

The U(N) symmetry can be reduced to SU(N) 
via the determinant only when one works with ex- 
actly N rishons on each link. This corresponds to 
working with the (2A^)!/(A r !) 2 -dimcnsional repre- 
sentation of SU(2N). Inserting the rishon repre- 
sentation of U x ,^ in the Hamilton operator, one 
can describe the quantum link dynamics as a hop- 
ping of rishons. This is illustrated in fig. 4. 
The partition function 



Z = Trexp(-/3H) 



(24) 
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of a quantum link model defined on a 4-d lattice 
can be expressed as a (4+l)-d path integral. Note 
that we have not imposed the 5-d Gauss law, i.e., 
we have not included a projector on gauge in- 
variant states propagating in the fifth direction. 
This implies that the fifth component of the non- 
Abelian vector potential, A$ — 0, vanishes. This 
is important, because it leaves us with the correct 
field content after dimensional reduction. Note 
that the physical 4-d Gauss law is properly im- 
posed, because the model does contain non-trivial 
Polyakov loops in the Euclidean time direction. 

For quantum spin models, dimensional reduc- 
tion arises naturally, because the spontaneous 
breakdown of a continuous global symmetry pro- 
vides massless Goldstone modes, and hence an 
infinite correlation length. On the other hand, 
when a gauge symmetry breaks spontaneously, 
due to the Higgs mechanism there are no massless 
modes, and hence dimensional reduction would 
not occur. Also in a confined phase the cor- 
relation length is finite. However, non-Abelian 
gauge theories in five dimensions generically have 
a massless Coulomb phase f2^| . Based on uni- 
versality, it is natural to assume the presence of 
a massless Coulomb phase also for quantum link 
models. Then the 5-d massless Coulombic gluons 
are described by the low-energy effective action 



S[A] = / dare 
Jo 



d*5 



X ^ TlF ^ F ^ 



-Tr dsApdsApl (25) 



which is just the standard 5-d Yang-Mills action 
with A$ = 0. The quantum link model is charac- 
terized by the "velocity of light" c, which is the 
analog of the spin- wave velocity in the Heisenberg 
model. Note that /u in eq. (|25|) runs over 4-d in- 
dices only. The dimensionful 5-d gauge coupling 
1/e 2 is the analog of p s in the spin model. At 
finite (3 the above theory has a 4-d gauge invari- 
ance only, because A^ = 0. At (3 — oo we are 
in the 5-d Coulomb phase with massless gluons, 
and hence with an infinite correlation length £. 
At finite (3 the extent of the extra dimension be- 
comes negligible compared to £, and the theory 
appears to be dimensionally reduced to four di- 
mensions. However, in four dimensions the con- 



finement hypothesis suggests that gluons are no 
longer massless. Then, as argued in refs.fflj^], a 
finite correlation length 



C K 6XP( ll^ 



(26) 



is generated non-perturbatively. Here 247r 2 /lliV 
is the 1-loop /^-function coefficient of SU(N) 
gauge theory, and 



1/5 2 = (3/e 



(27) 



is the gauge coupling of the dimensionally re- 
duced 4-d theory. The continuum limit g — ► 
of the 4-d theory is reached when the extent (3 
of the fifth direction becomes large. Like the 
spin model, one can consider the dimensionally 
reduced 4-d theory as a Wilsonian lattice theory 
with lattice spacing f3c (which has nothing to do 
with the microscopic lattice spacing of the quan- 
tum link model). As before, one averages the 5-d 
field over cubic blocks of size (3 in the fifth direc- 
tion and of size (3c in the four physical space-time 
directions. The block centers form a 4-d space- 
time lattice of spacing /3c, and the effective the- 
ory of the block averaged 5-d Coulombic gluons 
resembles a 4-d Wilsonian lattice gauge theory. 

Higher dimensions play an important role also 
in string theory. For example, the five known su- 
perstring theories are anomaly free only in ten di- 
mensions. The unifying low-energy effective the- 
ory of all these models is supergravity in eleven di- 
mensions. The various string theories are just dif- 
ferent 10-dimensional reductions of the same 11- 
dimensional effective theory. The coupling con- 
stant of the 10-d string theory is related to the 11- 
d supergravity coupling and to the extent of the 
extra dimension. Attempts to formulate string 
theory beyond perturbation theory have led to 
M-theory, which is a microscopic structure under- 
lying the non-renormalizable 11-dimensional su- 
pergravity. A candidate for M-theory is a matrix 
model [|||, in which the classical string coordi- 
nates are replaced by non-commuting oper- 
ators. All this is remarkably similar to the D- 
theory formulation of QCD in four dimensions. 
First, the classical gluon field A^ is replaced by 
non-commuting quantum link operators. The 
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low-energy effective theory of the resulting D- 
theory is a non-renormalizable 5-d non-Abelian 
gauge theory. Finally, just as in string theory, 
the resulting gauge coupling of QCD in four di- 
mensions 1/g 2 = f3/e 2 is related to the 5-d gauge 
coupling e and to the extent j3 of the fifth direc- 
tion. Due to these striking similarities, it seems 
that D-theory can be viewed as a daughter of M- 
theory — the string theorist's candidate for the 
mother of all fundamental theories. 

5. Numerical Verification of Dimensional 
Reduction in Wilsonian SU(3) Lattice 
Gauge Theory 4 

The above scenario of dimensional reduction 
crucially depends on the presence of a massless 
Coulomb phase in 5-d non-Abelian gauge theo- 
ries. Early numerical evidence for a weak cou- 
pling phase in 5-d SU(2) gauge theory was pre- 
sented in ref . . Here we extend these results to 
SU(3) and establish that the weak coupling phase 
is indeed a massless Coulomb phase 1 24 . The cor- 



responding numerical simulations are performed 
using Wilson's formulation of lattice gauge the- 
ory. Via universality they should also apply to 
D-theory. Once cluster algorithms become avail- 
able for non-Abelian quantum link models, the 
assumption of universality will be tested in de- 
tail. 

First, 5-d pure SU (3) lattice gauge theory has 
been simulated using the standard Wilson action 
on 4 5 — 8 5 lattices. There is a strong first order 
phase transition separating the confined phase at 
strong coupling from a weak coupling phase. The 
phase transition is at 6/e 2 , = 4.35(15). The ques- 
tion arises if the weak coupling phase is a mass- 
less Coulomb phase or a massive Higgs phase. To 
investigate this question, a Higgs model includ- 
ing an explicit scalar degree of freedom has been 
studied. In the unitary gauge, the action of the 
extended model is given by 



S[U] 



3 □ 

J2(l - ^Re TrU Xilt ). 



(28) 



Confined Phase 



0.5 

(6/e 2 ) / [ (6/e 2 ) + 4.0 ] 



Higgs Phase 



| | | | - 



Coulomb Phase 



Figure 5. The phase diagram of the 5-d SU(3) 
Higgs model: The weak coupling Coulomb phase 
is separated from the Higgs-confmcrncnt phase by 
strong first order phase transitions. 

The phase diagram of the model is shown in fig. 5. 
The weak coupling phase of the pure gauge theory 
extends to small values of k, but is separated from 
the Higgs phase at large k by a strong first order 
phase transition. The massive Higgs and confined 
phases are analytically connected. 

To establish that the weak coupling phase 
of the pure gauge theory is indeed a massless 
Coulomb phase, the static quark potential has 
been investigated on a 16 5 lattice at 6/e 2 , = 4.6 
and 9.0 for k = 0. Wilson loops W(R,T) of var- 
ious sizes (R, T) are measured and fitted to the 
functional form 



W(R,T) = Ae*p(-e 2 F(R,T)), 



(29) 



Based on a talk given by D. Chen 



where F(R,T) is the lattice tree-level perturba- 
tive result. Fig. 6 shows the fit for 6/e 2 , = 4.6. The 
fit parameters are A = 1.103(5), e 2 = 2.525(2), 
and x 2 /d.o.f. = 1.2(7). Similarly, for 6/e 2 , = 9.0 
one obtains A = 1.0089(3), e 2 = 0.8333(2), and 
X 2 /d.o.f. = 0.8(2). As it should, the parameter 
A approaches 1 in the eo — > limit. Also note 
that the charge renormalization e 2 /e^ drops from 
1.936 at 6/e 2 = 4.6 to 1.25 at 6/e 2 , = 9.0. Fi- 
nally, when a massterm for the gauge bosons is 
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Figure 6. Fit of Wilson loops: Loops of sizes 4x4 
to 8 x 8 on a 16 5 lattices at 6/eQ = 4.6 arc fitted 
to a 5-d lattice Coulomb potential. 

included in the fit, its fitted value is consistent 
with zero. This confirms the presence of a mass- 
less weak coupling phase with a 5-d Coulomb po- 
tential V(R) = e 2 /R 2 . 

To verify reduction from five to four dimen- 
sions, the 4-d finite temperature phase transition 
has been investigated using 5-d pure gauge the- 
ory with A$ = 0. The extent (3 in the fifth di- 
rection generates an effective 4-d gauge coupling 
l/g 2 — (3/e 2 . The extent in the fourth (Euclidean 
time) direction is also finite, and corresponds to 
the inverse physical temperature. For practical 
purposes, it is necessary to work with asymmet- 
ric lattices with a coupling 6/ep for the space- 
time plaquettes and a coupling 6/e 2 for plaque- 
ttes involving the fifth direction. The simulations 
were done on 16 3 x 6 x 2 lattices with 6/e§ = 4.6 
and 6/e§ = 2.41, 2.49. The Polyakov loop in 
the fourth (Euclidean time) direction is shown in 
fig. 7. One sees a clear signal for a confinement- 
deconfincmcnt phase transition in the dimension- 
ally reduced 4-d theory. Of course, when one uses 
Wilson's method, in the pure gauge theory there 
is no reason to go to five dimensions. Still, our 
results show that the expected dimensional reduc- 
tion to ordinary 4-d physics actually takes place, 



Figure 7. Behavior of the 4-d Polyakov loop: 
Varying the coupling &/e 2 of plaquettes involv- 
ing the fifth direction, a signal for a confinement- 
deconfinement phase transition is seen in the di- 
mensionally reduced 4-d theory. 

when one chooses the coupling such that the 5-d 
system is in the Coulomb phase. If the quantum 
link model falls in the Coulomb phase, univer- 
sality implies that it also undergoes dimensional 
reduction. Once cluster algorithms become avail- 
able for non-Abelian quantum link models, this 
scenario will be tested in detail. If powerful clus- 
ter algorithms can be constructed, it is likely that 
working with a 5-d quantum link model is more 
efficient than using Wilson's method in four di- 
mensions. 

6. Quantum Link QCD with Quarks 5 

Let us now construct full QCD in the frame- 
work of D-theory. The construction principle of 
D-theory is to replace the classical action in Wil- 
son's formulation by a Hamilton operator. For 
fcrmions one also replaces "4> x by 'J J. 75. Thus, the 
quantum link QCD Hamilton operator takes the 
form 

H = J ]T Tr{U x ,,U x+ ^Ul + ^UU 

5 Based on a talk given by R. C. Browcr 
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+ J'J2 [detC/^+detC/t M ] 

^ 1 + /i ^5 7m ^ a; ] 

a;, /j, 

+ M ^ *t 75 ^ + [2*t 75*, 

x a;, /j, 

- >l'> - *i +/> 7»C4, M *x]. (30) 

Here and ^ are quark creation and annihi- 
lation operators with canonical anticommutation 
relations 

{K aa ,% m } = 5 xy 6 ij S ab S al3 , 

v v »p} = {*r f , n m } = °> ( 3l ) 

where (a,b) and (a,/?) are color, flavor and 

Dirac indices, respectively. The generators of 
SU(N) gauge transformation now take the form 

G x = X)(^-a,m + + *£A*x. (32) 

The dimensional reduction of fermions is not 
completely straightforward. When one uses an- 
tiperiodic boundary conditions in the extra di- 
mension, the Matsubara modes, p§ = 27r(n5 + 
i)//3, lead to a short 0((3) fermionic correlation 
length. Hence, from the point of view of the 

4- d gluon dynamics, fermions with antiperiodic 
boundary conditions in the fifth direction stay at 
the cut-off and do not undergo dimensional re- 
duction. When one uses periodic boundary con- 
ditions, a Matsubara mode, = 0, arises, and 
the quarks survive dimensional reduction, but we 
face the same fine-tuning problem that arises for 
Wilson fermions. 

The fine-tuning problem has been solved very 
elegantly in Shamir's variant (25) of Kaplan's 
fermion proposal pq] . Kaplan coupled 5-d 
fermions to a 4-d domain wall, and found that 
a fermionic zero-mode gets bound to the wall. 
From the point of view of the dimensionally re- 
duced theory, the zero-mode represents a 4-d chi- 
ral fermion. For QCD, Shamir has simplified Ka- 
plan's proposal by formulating the theory in a 

5- d slab of finite size (5 with open boundary con- 
ditions for the fermions at the two sides. This ge- 
ometry fits naturally with the D-theory construc- 



tion of quantum link models. With open bound- 
ary conditions for the quarks and with periodic 
boundary conditions for the gluons, the partition 
function takes the form 

Z = Tr(0|cxp(-/3H)|0). (33) 

The trace extends over the gluonic Hilbert space 
only. Taking the expectation value in the Fock 
state |0), implies that there are no left-handed 
quarks at x§ = 0, and no right-handed quarks at 

In the presence of quarks, the low-energy effec- 
tive theory of the gluons (with A$ = 0) must be 
modified to 




x {^[TtF^F^ + -jTrdsApdsAp] 

+ ^(A„ + d IM )+M+^ 5 d 5 M. (34) 

The "velocity of light" c' of the quarks in the 
fifth direction is expected to be different from the 
velocity c of the gluons, because in D-theory there 
is no symmetry between the four physical space- 
time directions and the extra fifth direction. This 
is no problem, because we are only interested in 
the 4-d physics after dimensional reduction. 

Due to confinement, after dimensional reduc- 
tion the gluonic correlation length is exponen- 
tially large, but not infinite. As explained in 
ref.0, the same is true for the quarks, but 
for a different reason. Even free quarks pick 
up an exponentially small mass due to tunnel- 
ing between the two boundaries of the 5-d slab. 
The corresponding tunneling correlation length is 
2Mexp{M0). This suggests how D-theory can 
avoid the fine-tuning problem that arises for Wil- 
son fermions. In the chiral limit of quantum link 
QCD, the 4-d gluon dynamics takes place at a 
length scale 

24^ 2 /3 

^ eM (llN-2N f )e^ (35) 

which is determined by the 1-loop coefficient of 
the /3-function of QCD with Nf massless quarks 
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and by the 5-d gauge coupling e. If one chooses 
24tt 2 



M > 



(UN -2N f )e 2 



(36) 



the chiral limit is automatically approached be- 
fore one reaches the continuum limit as (3 becomes 
large. This solves the fermion doubling problem 
as well as the resulting fine-tuning problem in the 
D-theory formulation of QCD. 

7. A Flux Cluster Algorithm for the U(l) 
Quantum Link Model 6 

Cluster algorithms have been extremely suc- 
cessful in solving lattice field theories with a 
global symmetry, like Ising and Potts models |2q| 
as well as O(N) models |(]J2§] and $ 4 theory gCj. 
In these cases critical slowing down is practically 
eliminated, which allows very precise numerical 
simulations close to the continuum limit. In these 
cases, the clusters can be identified with physi- 
cal excitations |^l|p2[ . The cluster size is related 
to a susceptibility, which implies that the clus- 
ters cannot grow beyond the physical correlation 
length. This is essential for the efficiency of these 
algorithms. For other models, like CP(N) mod- 
els with N > 1 @ and SU(N) ® SU(N) chiral 
models with N > 2, cluster algorithms can also be 
constructed. However, they turn out not to be ef- 
ficient, as discussed in Ref. p3|. In these cases the 
clusters are not related to a physical quantity. For 
lattice gauge theories the situation is even more 
difficult. Only for models with discrete gauge 
groups efficient cluster algorithms could be found 
]35|j3(| . Despite numerous attempts, so far no 
efficient cluster algorithm has been constructed 
for models with continuous gauge groups. These 
attempts failed due to the presence of frustrated 
interactions and difficulties to find an efficient em- 
bedding of discrete variables. 

Here we construct the first efficient cluster algo- 
rithm for a model with a continuous gauge group 
- the U(l) quantum link model. When for- 
mulated as a (d + l)-dimensional D-theory, this 
model is expected to reproduce the physics of 
the standard Wilson formulation of e?-dimensional 
compact E/(l) gauge theory. For d = 4 there 

6 Based on a talk given by A. Tsapalis 



is a phase transition between a strong coupling 
confined phase and a weak coupling massless 
Coulomb phase |p7| . Due to the infinite correla- 
tion length of 5-d photons, dimensional reduction 
from 5 to 4 dimensions occurs when the extent (3 
of the additional fifth dimension becomes finite. 
As in non-Abelian quantum link models, one has 
1/g 2 = (3/e 2 . Dimensional reduction occurs as 
long as the effective 4-d gauge coupling g is weak 
enough to fall in the Coulomb phase of the 4- 
d theory. At stronger coupling, and hence at a 
smaller extent of the fifth direction, one enters 
the confined phase. Then the correlation length 
is finite, and dimensional reduction does not oc- 
cur. For d = 3 compact U(l) gauge theory al- 
ways confines |38|,|39|], but the correlation length 
diverges exponentially in the weak coupling limit. 
Then the mechanism of dimensional reduction is 
very similar to the one in QCD going from 5 to 4 
dimensions. In particular, dimensional reduction 
occurs only when the extent of the extra dimen- 
sion becomes large. 

An efficient cluster algorithm for the U(l) 
quantum link model can be constructed, be- 
cause the variables are discrete even though the 
gauge symmetry is continuous. Hence, in contrast 
to Wilson's formulation, no embedding problem 
arises. Also, the clusters are physical objects - 
namely world-sheets of electric flux strings prop- 
agating in the additional Euclidean dimension. 
The flux cluster algorithm is a gauge analog of 
the loop algorithm for the quantum Heisenberg 
model |jl^,0. Like the loop algorithm, it oper- 
ates directly in the continuum of the additional 
Euclidean dimension || . Improved estimators for 
Wilson loops can also be constructed. 

The Hamilton operator of the U(l) quantum 
link model is defined on a d-dimensional lattice 
and is given by 

2 ^ 

X,fJb<U 



[u x ,,u x+ ^ul +0 M 



+ U X M X 



,U1 



.ul 



(37) 



The quantum link operators U Xjfl and the gener- 
ators of gauge transformations G x can be repre- 
sented as 

U x ,t_i = S xfl + iS X)i — S x>/M , 
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where S 1 * „ is a spin operator associated with a 
link, with the usual commutation relations 

= i$x,vf>iui£ijkS Xtli . (39) 

The operator = T x ^/2 represents the elec- 
tric flux, and L^.^ acts as a flux raising operator. 
The partition function 

Z = Trexp(-/lff) (40) 

can be formulated as a (d + l)-dimensional path 
integral. Here we describe the situation in d = 2. 
Performing a checker board decomposition H = 
H e + H a of the Hamiltonian into contributions 
of even and odd plaquettes, and discretizing the 
extra dimension into N intervals of size e, the 
partition function takes the form 

Z = lim Tr{[exp(-e# e ) exp(~eH )] N }. (41) 

Inserting complete sets of basis states between 
each pair of operators yields a (2 + l)-dimensional 
lattice with 2N slices in the third dimension. 
Here we work in an electric flux basis, and we 
limit ourselves to spin 1/2, such that the eigenval- 
ues of S XjfJi in the slice n (with n £ {1, 2, 2iV}) 
at t = en/2 are limited to e Xifli t = ±1/2. The 
partition function then takes the form 

Z =Ii E cxp(-S[e]). (42) 

The action S[e] is a sum of contributions from 
checker boarded cubes carrying an 8-link interac- 
tion. The Boltzmann factor of a cube is deter- 
mined by the elements of the 16 x 16 plaquette 
transfer matrix 

r = cx P(e 2 l U ^ U ^+^ U l+u^ U x,u 

+ U*,»U x+ ^Ul + ^Ulj). (43) 

All diagonal elements of the transfer matrix are 

1, except (+ H |T| + + --) = ( \- + \T\ - 

— h +) = cosh(eJ/2) = c, and all off-diagonal 

elements vanish, except (+ H \T\ h +) = 

(- - + + |T| + + - -) = sinh(e J/2) = s. Here, ± 
denotes the electric fluxes e x ^, e x +p, jt/ , e x +p^, 



e x ,v around a plaquette in two adjacent time- 
slices. 

The cluster algorithm constructs open or closed 
oriented surfaces, which represent world-sheets of 
electric flux strings propagating in the extra di- 
mension. An update reverses the flux of all links 
belonging to the cluster. Cluster growth starts 
by selecting an initial link (x,/j,,t) at random. 
The flux variable e Xtflt t participates in two cube 
interactions, one at earlier, the other at later t 
values. The eight links of a cube are connected 
to clusters according to the following rules. For 
a cube configuration of weight c, with a proba- 
bility p — (c — s + 7) /8c each link is connected 
to its i-partner shifted in the extra dimension. 
With probability 1 — p the links are connected 
in two groups of four links, such that the four 
links with the same t value are connected. For 
a configuration of weight 1, with a probability 
q = (c — s + 7)/8 each link is connected to its 
i-partner. With probability 1 — q the links are 
connected in two groups of four, such that after 
flipping one group, a configuration of weight s is 
obtained. For a configuration of weight s, with 
probability r = (7c + s — 7) /8s the four links with 
the same i-value are connected. With probabil- 
ity (1 — r)/7 they are connected in two groups of 
four, such that after flipping one group, one of 
the fourteen weight 1 configurations is obtained. 
The generated clusters form open or closed sur- 
faces. The open surfaces arise, because we have 
not imposed the 5-d Gauss law for states prop- 
agating in the extra dimension. The algorithm 
is constructed such that it obeys ergodicity and 
detailed balance. 

Due to the discrete basis of the Hilbert space, 
the path integral expression for the partition 
function of eq.(|4l|) is well-defined in the con- 
tinuum of the extra dimension, i.e., at e = 0. 
Taking the continuum limit of the cluster rules, 
one obtains an algorithm that operates directly 
in the continuum of the extra dimension. This 
algorithm has been implemented along the lines 
of the continuous-time loop-cluster algorithm for 
the quantum Heisenberg model § . 

To show that the clusters are directly related 
to physical objects, we consider the expectation 
values of Wilson loops Wc — the product of quan- 
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turn link operators U Xtfl along some closed curve 
C on the d-dimensional lattice at a fixed value of 
t. The quantum links act as raising operators of 
electric flux. Hence, a configuration contribut- 
ing to (We) is inconsistent with the constraints 
in the configurations contributing to Z. Thus, it 
seems that just using the flux cluster algorithm 
does not provide information about Wilson loops. 
However, one can imagine flipping only a part of 
a cluster, which is bounded by a curve C. This 
leads into the right sector for collecting informa- 
tion about (Wc). In practice, one need not even 
perform these partial cluster flips explicitly, be- 
cause one can extract the same information by 
examining the clusters of the original algorithm. 
This leads to an improved estimator for Wilson 
loops, which only receives positive contributions 
from each cluster, and shows that the clusters 
themselves represent physical objects. To verify 
the efficiency of the algorithm we have measured 
autocorrelation times, which indeed are small. A 
detailed study of the dynamical critical exponent 
z is in progress. 

The flux cluster algorithm operating in the 
continuum of the extra dimension has been ap- 
plied to the (4 + l)-dimcnsional and (3 + 1)- 
dimcnsional C/(l) quantum link models. In the 
(4 + l)-dimensional model the cluster size per 5- 
d volume |C|/V is shown as a function of the 
extent (3 of the fifth direction in fig. 8. Above 
Pc sa 0.65 the cluster size is proportional to the 
volume, indicating the infinite correlation length 
in the Coulomb phase. As shown in fig. 9, in the 
(3 + l)-dimensional model, |C|/V goes to zero 
when one increases the volume, which is consis- 
tent with a confined phase for all f3. As expected, 
a large increase of the cluster size is observed in 
the large j3 limit. A detailed investigation of Wil- 
son loops in d = 3 and d = 4 using improved 
estimators is in progress. 

8. Conclusions 

D-theory provides a new non-perturbative for- 
mulation of quantum field theory. Dimensional 
reduction of discrete variables is a generic phe- 
nomenon that occurs for various models, includ- 
ing O(N), CP(N), and SU(N) ® SU(N) scalar 
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Figure 8. Cluster area per 5-d volume in the 
(4 + l)-d U(l) quantum link model: Above the 
critical extent (3 C ~ 0.65 of the fifth direction, the 
cluster size increases with the volume, indicating 
the presence of a Coulomb phase with an infinite 
correlation length. 

models, as well as Abelian and non-Abelian gauge 
theories, in particular QCD. The discrete nature 
of the fundamental variables makes D-theory at- 
tractive, both from an analytic and from a com- 
putational point of view. On the analytic side, 
the discrete variables allow us to rewrite the 
bosonic fields in terms of fermionic rishon con- 
stituents. This may turn out to be useful for 
studying the large N limit. On the numerical 
side powerful cluster algorithms become available, 
which may dramatically improve numerical sim- 
ulations of lattice field theories. A lot of work 
needs to be done to decide if D-theory provides 
a more efficient non-perturbative quantization of 
field theories than Wilson's method. 
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